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SUMMARY 


The NASTRAN preprocessor BANDIT, which improves NASTRAN's computer 
efficiency by resequencing grid point labels for reduced matrix bandwidth, 
has been improved by the addition of (1) the Gibbs-Poole-Stockmeyer (GPS) 

I algorithm, and (2) the user option to reduce matrix profile rather than 
I matrix bandwidth. After describing these program additions, this paper 
I shows that, compared to the Cuthill-McKee (CM) algorithm on which BANDIT 
was originally based, GPS is faster and achieves similar results. For 
I completeness, BANDIT’ s current capabilities and options are summarized. 


BACKGROUND 


The NASTRAN structural analysis computer program (ref. 1, 2), as a 
finite element program, assembles matrices which are normally both symmetric 
and sparsely-populated. The locations of the nonzero terms in the matrices 
are determined solely by the choice of numbers (labels) assigned to the grid 
I points. Like most finite element codes, NASTRAN’s computer running time 
can be reduced if the labels can be chosen in such a way that the nonzero 
terms cluster tightly about the main diagonal. The NASTRAN user has com- 
plete control over that clustering by hia choice of grid point labels and 
I his optional use of SEQGP bulk data cards, which effect an internal grid 
point resequencing for calculation purposes. 

Soon after NASTRAN became available some five years ago, it was 
I apparent that the program user could benefit from an automatic capability 
: to perform the resequencing and generate the SEQGP cards. Indeed, for 
large complex structures or those generated automatically, the job of 
determining a good grid point sequence manually was, at best, tedious and 
often very difficult. 

To fill the need for an automatic capability, several NASTRAN pre- 
' processor computer programs were developed: BANDIT (refs. 3, 4), WAVEFRONT 
(refs. 5-7), and BANDAID (ref. 8). (For a general survey of NASTRAN pre- 
processors and postprocessors, see reference 9.) Both BANDIT and BANDAID 
are intended to reduce matrix bandwidth, while WAVEFRONT is intended to 
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reduce matrix wavefront. (These terms are defined in the next section.) Of 
these preprocessors, BANDIT and WAVEFRONT appear to be the most popular. 
BANDIT was originally based on the Cuthill-McKee resequencing algorithm 
(ref. 10). WAVEFRONT and BANDAID are based on strategies developed by their 
authors. Levy (ref. 5) and Cook (ref. 8), respectively. These algorithms 
and others have been reviewed and compared by Cuthill (ref . 11) . 

Recently, a new bandwidth and profile reducing algorithm was developed 
by Gibbs, Poole, and Stockmeyer (GPS) (ref. 12) of The College of William 
and Mary. Since their testing of it showed it to be both effective and 
efficient (ref. 13), we have incorporated it in the BANDIT program to 
supplement the Cuthill-McKee (CM) strategy already there. (Actually, BANDIT 
uses the so-called reverse Cuthill-McKee algorithm since it was observed 
by George (ref. 14) and later proved by Liu and Sherman (ref. 15) that 
reversing the sequence generated by CM can never increase the profile and 
frequently reduces it. Such a reversal has no effect on the matrix band- 
width.) In general, GPS executes faster than CM and achieves comparable 
results. Unfortunately, for a given structure, it is not possible to predict 
a priori which strategy will yield the smaller matrix bandwidth or profile. 
However, since excessive resequencing time has never been considered to be a 
problem, BANDIT 's current default mode of operation is to apply both CM and 
GPS to the structure in order to get the better of the two results . 


DEFINITIONS 


For the purposes of this discussion, some useful terms will be defined 
which generally follow the material given in Cuthill 's survey (ref. 11). 

Given a symmetric matrix A of order N, we define a "row bandwidth" b^ 
for row i to be the number of columns separatxng the first nonzero in the 
row from the diagonal. Alternatively, b^^ is the dxfference between i and 
the column index of the first nonzero entry of row i of A. Then the matrix 
bandwidth B and profile P are defined as 

B = max b , (1) 

i<N ^ 

N 


Let w. denote the number of active columns in row i. A column j is 
active in row i if j>i and there is a nonzero entry in that column in any 
row with index kii. Thus, a given column is activated at the first nonzero 
encountered (reading from top to bottom) and remains active until the 
diagonal is reached. The matrix wavefront W is then defined as 


(3) 
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W = max w. 
i<N ^ 


(4) 


I 

i 

I 

Since the matrix A is symmetric, 

N N 

j P = E b. = E w. 

' i=l ^ i=l ^ 

I Now, for row 1, let ^ denote the columnar distance between the diagonal 
and the last active column in row i. Then 

! 

' B = max b. = max b, (5) 

' i^N ^ ii^N 

Since, by definition. 


j for each i, it 


and 


Hence, as a consequence of these definitions, the matrix wavefront W 
^ for a given matrix is less than the matrix bandwidth B, and the matrix 
I profile P is equal to both the sum of the "row bandwidths" and the sum of 
the "row wavefronts." 

i 

These definitions are generally modified slightly by preprocessors such 
as BANDIT. Since NASTRAN requires all external resequencing via SEQGP cards 
to be performed at the grid point level rather than the degree of freedom 
(DOF) level, BANDIT treats each grid point as if it had only one DOF. In 
general, a NASTRAN grid point can have as many as six DOF’s. Thus, to 
convert BANDIT' s values of bandwidth and profile to meaningful approximate 
values for NASTRAN' s structural matrices, one must multiply by the average 
number of DOF's per grid point. 


b ~ w 


( 6 ) 


^s 


that 

B 


max h, - max w. = W 
i^N ^ i^N ^ 


(7) 


N N 

S = E b, - E w, = P 
1=1 1=1 ^ 


( 8 ) 


A NEW RESEQUENCING STRATEGY 


The principal recent improvement to BANDIT is the installation of the 
new bandwidth and profile reducing algorithm developed by Gibbs, Poole, and 
Stockmeyer (GPS) (refs. 12, 16) of The College of William and Mary. Rather 
than describe how GPS works, we shall instead demonstrate its performance 
on a set of test problems. The test problems used here constitute the 
current extent of a growing collection of diversified NASTRAN data decks to be 
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used for the testing of resequencing and equation solving algorithms. It is 
expected that a complete description of the set, including plots of each 
structure, will eventually be published. 

The results of the resequencing tests are shown in Table 1. In that 
table, the following definitions apply: 


N 

M 

B 

P 

T 

Orig. 

CM 

GPS 

Decomp . 


number of grid points (nodes) 

mavImuTTi nodal degree (i.e., the maximum number of nodes 
connected to any node) 

matrix bandwidth (in terms of grid points rather than DOF) 
matrix profile (in terms of grid points) 
time, CDC 6400 CP seconds 

an original value (before resequencing) of B or P 
Cuthill-McKee strategy 
Gibbs-Poole-Stockmeyer algorithm 
matrix decomposition 


For each of 20 structures, ranging in size up to 2680 grid points, the 
grid point labels were resequenced using both CM and GPS. Before and after 
results for both bandwidth (B) and profile (P) are shown. Since the test 
criterion was to reduce B rather than P, the P results are less significant. 
With CM, a user choice of profile reduction rather than bandwidth reduction 
will generally give different results for both P and B. All tests were run 
on a CDC 6400 computer with the SCOPE 3.4.2 operating system. Central 
processor (CP) times are given for both CM and GPS. 

Since some of the structures are clearly very large, rough estimates of 
the NASTRAN real, symmetric, single-precision decomposition times on a 
CDC 6400 are given in the last column of Table 1. These values were 
computed using the following formula extracted from the NASTRAN subroutine 
RSPSDC : 


T = Tg(n-2b/3)b^/2 + Tp(n-b/2)b (9) 

For decomposition times in Table 1, it is assumed that (1) there are no 
active columns (in the NASTRAN sense), (2) no "spill" occurs, and (3) the 
structure has six DOF’s per node. Hence n=6N and b=6B, where the bandwidth 
B used is the minimum of that obtained by CM and GPS. The constants T^ and 
Tp are computer-dependent time constants equal, respectively, to 15 ysec and 

140 ysec for the CDC 6400. 
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Several conclusions can be drawn from the table: 

1. CM and GPS generally obtain comparable bandwidth results, 
although occasionally one does significantly better than the other. 

2. GPS is faster than CM. 

3. Both CM and GPS are generally fast compared to estimated 
decomposition times. In the absence of resequencing, the decomposition 
times would usually be much larger. 

Conclusions 1 and 3 indicate that the user would, in general, benefit 
from having both CM and GPS attempt to resequence his structure. Thus, the 
default mode of operation in BANDIT uses both and delivers to the user SEQGP 
cards for the better result. 


REDUCTION OF MATRIX PROFILE 


The second recent improvement to BANDIT is that the user now has the 
option of selecting matrix profile reduction rather than matrix bandwidth 
reduction. This option was installed primarily to facilitate testing with 
NASTRAN Level 15.9 to determine whether profile reduction has any 
advantages over band reduction. At this writing the question is still open. 
However, based on the close relationship between a matrix's bandwidth and 
its profile, it seems unlikely that major advantages wxll result. Indeed, 
in a larger sense, equation solvers which exploit matrix bandwidth, profile, 
or wavefront can all be classified under the general category of "envelope 
methods" (ref. 15), which ignore only those zeros in a matrix outside a 
particular region of the matrix. Distinct from the envelope methods are the 
general sparse methods, which ignore all the zeros in a matrix. 


CURRENT BANDIT USAGE 


This section summarizes briefly how a NASTRAN user runs BANDIT and what 
bandit's list of options are. It is assumed here that the prospective 
BANDIT user has already compiled the program and has it in executable form. 

Versions of BANDIT exist for all computers on which NASTRAN runs: 

CDC 6000, IBM 360/370, UNIVAC 1100, and Honeywell 6000 (ref. 17). 

Input to BANDIT generally consists of a standard NASTRAN data deck (ID 
through ENDDATA) plus one or more special $ cards (which are comments to 
NASTRAN) for supplying various instructions to BANDIT. The minimum BANDIT 
data deck consists of $ option cards, BEGIN BULK, element connection cards, 
and ENDDATA. BANDIT does not use GRID cards. 
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Output from BANDIT consists generally of printed output, punched 
output, and a file (FORTRAN logical unit 8) containing the complete input 
I deck plus any SEQGP cards generated. This file, which is created 
automatically, is rewound before BANDIT execution terminates so that it is 
^ ready to be used as input to NASTRAN. 

I 

The current version of BANDIT, designated Version 5.1 and dated 
' 04/28/75, contains in its element library all NASTRAN elements in Level 15.5 
plus some additional elements appearing in several non-standard versions of 
NASTRAN. Multipoint constraint (MFC) cards are also recognized and 
I accounted for if the user so elects. 

' Instructions from the user to BANDIT are passed via $ cards having the 

general format 

^ $KEYW0RD1 KEYW0RD2 

i 

! where the $ must appear in card column 1, and the first letter of KEYWORDl 
j must appear in coltimn 2. Otherwise, the format of such cards is free field: 

' keywords, which can contain no embedded blanks, must be separated by one 
; or more blanks, and at least two letters of each keyword are required for 
recognition by BANDIT. Since the $ cards are interpreted by NASTRAN as 
I comments, they can be left in the deck during a NASTRAN run. 

j The complete list of current $ cards is summarized in Table 2. Such 

j cards can appear in any order but must be placed somewhere ahead of BEGIN 
! BULK. The cards defined under Part B are specialized cards created for 
' particular users with special needs. For most $ cards, a default is defined 
and denoted in Table 2 by underlining. The default applies whenever the 
$ card is omitted from the deck. 

I For example, referring to Table 2, if resequencing is to be performed, 

the user inserts the card 

$ SEQUENCE YES 

into the deck anywhere before the BEGIN BULK card. In most cases, this is 
' the only $ card added to the deck. 

! Although many of the cards listed in Table 2 are probably self- 

explanatory, several require additional explanation. The $GRID card is 
used to declare an upper bound (preferably least upper bound) on the number 
of grid points. The inclusion of this card is sometimes necessary (and 
never hurts) if BANDIT' s default allocation of "open core" to various tables 
is inadequate. Generally, the default is such that the maximum nodal degree 
I is limited to about 19. (The degree of a node is the number of other nodes 
j connected to it.) Thus, for example, a $GRID card is required whenever 
! solid elements are present. 

Sometimes, in order to induce active columns in NASTRAN, the user 
would like BANDIT to ignore connections to selected grid points. Such 



TABLE 2 - SUMMARY OF BANDIT $ CARDS 


A. For General Use 


$ SEQUENCE (^, YES) 

$PUNCH (NONE, SEQGP , ALL) 
$CRITERION ( BAND . PROFILE) 
$METHOD (CM, GPS, BOTH ) 
$MPC (NO, YES) 

$PRINT ( MIN . MAX) 

$GRID N 
$IGNORE G1,G2,--- 


Is resequencing to be performed? 
What should be punched? 

What should be reduced? 

By what method? 

Take MFC's into account? 

What printed output? 

Upper bound on number of grids. 
Grid points to ignore. 


B. For Particular Users 


$NASTRAN (^, YES) 

$ INSERT 
$ INSERT N 
$LINES N 
$PLUS + 

$CONNECTION (NO, YES) 
$START G1,G2,-*- 
$DEGREE N 


NASTRAN to follow BANDIT? 

Location of cards to insert. 

Number and location of cards to insert 
Number of lines per page. 

User-defined plus sign. 

Punch connection table? 

User-supplied CM starting nodes. 

Ignore nodes of degree exceeding N. 


! 


points are listed on the $ IGNORE card and are resequenced last. 

J The $MPC card is used to tell BANDIT to modify the matrix connectivity 

according to the multipoint constraints (MFC's) in the deck. If this option 
I is invoked, all MFC's present are included, regardless of any set ID's. 

The presence of MFC's creates a dilemma from the resequencing point— of— view, 

I since resequencing is always performed at the grid point level, whereas 
I MFC's are always applied at the DOF level. BANDIT treats MFC's by first 
, generating additional connections between each independent point in the 
constraint relation and every other point to which the dependent point was 
previously connected. Second, each dependent point is eliminated from the 
connection table. Thus, if most or all of the DOF's for the dependent 
points appear in MFC relations (as, for example, with rigid links), MFC's 
I should be taken into account. This guideline is based on experience with 
NASTRAN Level 15.5 and will probably have to be modified with Level 15.9 and 
; subsequent versions, since a new equation solver has been developed for 
! them (ref. 18). 

The $NASTRAN card was created for IBM users wanting to create a single 
j BANDIT— NASTRAN cataloged procedure in which the user could execute either 
j BANDIT or NASTRAN (or both) and be able to control the choice with $ cards . 
j The YES choice results in a FORTRAN STOF 5 at successful termination, thus 
supplying a testable condition code to the cataloged procedure. 


CONCLUDING REMARKS 


j From the test results presented, it is clear that the addition of the 

new resequencing strategy by Gibbs, Poole, and Stockmeyer greatly enhances 
I bandit's capabilities. The addition of the user option to reduce matrix 
profile rather than matrix bandwidth is a useful addition, but testing with 
NASTRAN Level 15.9 will be required to determine the extent of its usefulness. 
From the NASTRAN user’s point of view, the relevant question is: For 

Levels 15.9 and 16, how should the grid point labels be sequenced? When 
these versions become available, this question will hopefully be answered by 
testing with band, profile, and wavefront reducers. 
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